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INTRODUCTION 

A  thick-walled  cylinder  is  used  for  a  variety  of  applications  in  the 
chemical,  nuclear,  and  armament  industries  where  large  internal  pressures  have 
to  be  withstood.  In  the  absence  of  residual  stresses,  the  cracks  usually  form 
at  the  bore  where  the  hoop  stress  developed  by  the  pressure  is  highest.  To 
prevent  such  failure  and  to  increase  the  pressure-carrying  capacity,  a  common 
practice  is  autofrettage  treatment  of  the  cylinder  prior  to  use  (ref  1).  This 
process  has  the  effect  of  producing  beneficial  compressive  residual  hoop 
stresses  near  the  bore  which  can  prolong  the  fatigue  life. 

The  determination  of  residual  stresses  in  autof rettaged  cylinders  has  been 
considered  by  many  investigators  using  different  mathematical  methods  and 
material  models  (refs  2-8).  Most  of  the  earlier  solutions  were  based  on  the 
assumption  that  the  material  behaves  elastically  on  the  release  of  the 
autofrettage  pressure.  However,  many  materials,  particularly  the  quenched  and 
tempered,  low  alloy  steels  generally  used  for  high  pressure  vessels,  exhibit  a 
significant  Bauschinger  effect  (ref  9).  In  a  recent  paper  (ref  10),  this 
author  presented  a  closed-form  solution  of  residual  stresses  in  auto* rettaged 
tubes  based  on  a  theoretical  model  considering  the  Bauschinger  and  hardening 
effects  during  unloading,  but  neglecting  the  strain-hardening  effects  during 
loading.  A  more  general  theoretical  model  without  this  restriction  was  pro¬ 
posed  earlier,  but  only  part  of  the  final  results  were  shown  due  to  space 
limitation  (ref  11).  In  the  present  report,  the  complete  method  of  stress  and 
deformation  analysis  based  on  the  general  theoretical  model  is  stated.  The 
new  model  is  a  better  representation  of  the  actual  loading/unloading  behavior 

References  are  listed  at  the  end  of  this  report. 


in  a  high  strength  steel.  The  Bauschinger  effect  factor  is  treated  as  a  func¬ 
tion  of  overstrain.  The  strain-hardening  effect  is  taken  into  account  with 
different  parameters  used  for  loading  and  unloading  processes.  The  formulas 
for  calculating  stresses,  strains,  and  displacements  are  given  and  new  results 
of  residual  stresses  in  autofrettaged  high  pressure  vessels  are  presented. 

MATERIAL  BEHAVIOR  AND  MODELING 

The  material  chosen  for  this  investigation  was  a  modified  4330  steel 
having  a  martensitic  structure.  A  description  of  its  chemical  composition  and 
various  heat  treatments  is  given  in  Reference  9  by  Milligan,  Koo,  and 
Davidson.  They  studied  the  material  behavior  by  utilizing  a  uniaxial  tension- 
compression  specimen.  Figure  1  shows  the  stress-strain  curve  during  loading 
and  unloading  after  overstrains  in  tension.  The  stress-strain  curve  during 
loading  can  be  replaced  with  sufficient  accuracy  by  a  bilinear  elastic-plastic 
model  shown  in  the  figure  by  the  broken  lines.  For  the  plastic  portion,  the 
yield  stress  (a)  is  related  to  the  plastic  strain  (eP)  by 

a/aQ  =  l  +  mC/(l-m)  and  f  =  (E/o0)eP  (1) 

where  Young's  modulus  (E),  tangent  modulus  (mE),  initial  yield  stress  (aQ) , 
and  the  Poisson's  ratio  (v)  are  the  material  constants.  The  nominal  yield 
strength  at  0.1  percent  offset  was  chosen  as  the  initial  yield  stress. 

Initially,  the  yield  stresses  in  tension  and  compression  are  approximately 
equal  so  that  the  material  can  be  considered  as  isotropic.  However,  the  ratio 
of  the  yield  stress  upon  reverse  yielding  to  the  initial  yield  stress  is 
strongly  affected  by  overstrain  in  a  high  strength  steel.  The  values  of  the 
Bauschinger  effect  factor  (BEF)  also  depend  on  the  0.1  percent  offset  men¬ 
tioned  above.  Taking  into  account  strain-hardening,  the  definition  of  the  BEF 


(2) 


BEF  =  (a0'-al)/al  =  f(eiP) 

where  aj,  ejP  are  the  yield  stress,  plastic  strain  just  before  unloading 
occurs,  and  a0'  is  the  linear  drop  in  stress  until  reverse  yielding  begins. 
Figure  2  shows  the  Bauschinger  effect  factor  (f)  as  a  function  of  percent  ten¬ 
sile  overstrain  (eP).  The  graph  shows  a  decrease  of  the  BEF  with  increasing 
amount  of  tensile  prestrain  up  to  approximately  two  percent  at  which  point  it 
becomes  effectively  constant. 

In  most  of  the  plasticity  theories,  the  curve  of  reverse  loading  is 
uniquely  defined  by  the  curve  of  the  first  loading.  The  present  model  does 
not  assume  such  a  relationship.  The  experimental  stress-strain  curve  during 
unloading  is  used  directly.  A  piecewise  linear  representation  can  be  used, 
but  only  a  bilinear  approximation  was  chosen  here  as  shown  in  Figure  1.  In 
the  same  figure,  we  also  show  two  other  theoretical  models  in  dotted  lines  for 
this  material  subjected  to  unloading  with  further  plastic  flow.  They  are 
isotropic  and  kinematic  hardening  models.  A  summary  of  the  plasticity 
theories  given  by  Armen  (ref  12)  illustrates  the  fact  that  all  of  the  theories 
are  capable  of  treating  the  monotonic  loading  situation.  For  cyclic  loading, 
including  reverse  yielding,  kinematic  and  isotropic  hardening  represent  the 
limits  to  the  actual  behavior,  whereas  the  remaining  theories  are  falling 
anywhere  within  these  limits.  The  above  conclusion  indicates  that  none  of  the 
existing  models  can  accurately  represent  the  actual  material  behavior  as  shown 
in  Figure  1.  The  reverse  yielding  begins  at  a  stress  whose  magnitude  is 
smaller  than  that  predicted  by  all  existing  theories  and  the  slope  of  strain¬ 
hardening  after  reverse  yielding  is  much  larger  than  that  during  the  first 
loading  process.  The  new  bilinear  model  presented  here  is  a  better  approxima¬ 
tion  to  the  actual  material  behavior. 
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Figure  1.  Stress-strain  curve  during  loading  and  unloading. 
Choosing  a  new  coordinate  system  (ct’.e’)  with  the  origin  at  the  point 
before  unloading,  we  have  for  the  plastic  portion  of  the  reverse  yielding 


curve 


/°o  s  *o'/a0  +  m’C’/d-n*’ )  and  C'  =  (E/o0)e’P  (3) 

where  m'E  is  the  slope  of  the  reverse  yielding  curve  and  e'P  is  the  additional 
plastic  strain  during  unloading.  According  to  Eqs.  (1)  and  (2),  cr0 '  can  be 
expressed  as  a  function  of  plastic  strain  (EjP)  just  prior  to  unloading  by 

<V/<70  =  [1  ♦  mCi/(l-m)][l  +  f  (Ci )  ]  =  g  ( C  i )  (4) 

With  experimental  data  f(cP)  shown  in  Figure  2,  we  can  calculate  o0'  and 
determine  the  initiation  of  reverse  yielding.  The  slope  of  reverse  yielding 


subject  to  ag  >  az  >  <xr 


<?r/ao  1  p*  p2  (^) 

=-(--+  --)  for  p  <  r  <  b 

<70/<7o  2  b*  r*  (8) 

oz/o0  =  v(Pr+cr0)/cro  ♦  ^z/^o  (9) 

E  u  ar  p*  E 

-  =  (l-2v) (l+v)  —  +  (1-u2)  -  w  —  ez  (10) 

<*o  r  ao  *o 


and 

(M-2u) 

(E/o0)£2  .  -----  (P/<70)  (11) 


where  u  =  0  (open-end),  1  (closed-end),  and  p  is  the  elastic-plastic  boundary 
relating  to  the  internal  pressure  p  by 


1  1 

p/o0  =  -  (l-p*/b* )  ♦  (I-82) log(p/a)  +  -  02(p2/a2-l)  (12) 

2  2 


The  equivalent  plastic  strain  can  be  calculated  by 

(E/<70)cP  *  (  *  8i(P*/r*-l)  ,  in  (a  <  r  <  p)  (13) 

and 

1-m 

01  = -  -  ,  02  *  m8i/(l-m)  (14) 

ft  d-.i 

m  ♦  - 

2  ( 1-y*  ) 


REVERSE  YIELOING 

If  the  pressure  p  given  by  Eq.  (12)  is  subsequently  removed  completely 
with  no  reverse  yielding,  the  unloading  is  entirely  elastic  and  the  solution 
is  given  by 


°r' 

p  b* 

*  -----  [t  -  -  1) 

(15) 

<V 

b*/a*-l  r* 

(16) 

^z' 

■  y  (<7r  '  +<70  ' )  +■  E  £Z' 

(17) 

Eez' 

=  -  (p-2v)p/(b*/a*-l ) 

(18) 
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Eu'/r  =  -[(l-v-pp)  +  (1+y  )bVra]p/(b2/aJ-l)  (19) 

Let  a  double  prime  denote  a  component  in  the  residual  state,  i.e.,  og"  = 
oq  oq'  .  Assuming  a  reduced  compressive  yield  strength  as  a  result  of  the 
Bauschinger  effect  and  using  Tresca's  yield  criterion  subject  to  or"  >  az"  > 
oq"  ,  the  reverse  yielding  will  not  occur  if 

<7r"  -  a g"  *  o"  =  fa0[l  +  mC/(l-m)]  (20) 

Substituting  the  loading  and  unloading  solutions  into  Eq.  (20),  we  can  deter¬ 
mine  the  minimum  pressure  (pm)  for  reverse  yielding  to  occur.  The  equation 
for  pm  is  given  by 

1 

Pm/<7o  =  “  (l-a*/b«)[l  +  mCm/(l-m)][l  +  f (Cm) ]  (21) 

and  Cm  =  . 

Equating  Eq.  (12)  to  (21),  we  can  determine  the  minimum  amount  of 
overstrain  (pm)  and  calculate  the  minimum  pressure  (pm)  required  for  reverse 
yielding . 

ELASTIC-PLASTIC  UNLOADING 

Now  suppose  that  the  loading  has  been  such  that  the  internal  pressure  is 
larger  than  pm  given  by  Eq.  (21).  On  unloading,  yielding  will  occur  for  a  <  r 
<  p'  with  p'  <  p.  In  the  reverse  yielding  zone,  the  stresses  in  the  residual 
state  satisfy 

<7r"  -  oq"  =  o"  (22) 

assuming  that  crr"  J  oz"  >  erg".  Since  the  yield  criterion  during  loading  is 

o9  -  ar  =  a  (23) 

we  have 

CTp'  -  Oq'  >  (<7p"-C7g")  +  (<70-ar)  =  o"  +  O  =  CT ’  (24) 

where  o  and  o'  are  defined  by  Eqs.  (1)  and  (3),  respectively. 


The  associated  flow  rule  states  that 


de0'P  =  -der'P  $  0  and  de2’P  =  0  (25) 

Therefore,  we  have 

ez'  =  ez'P  =  -£z  (26) 

Since  ez'  is  now  known,  Hooke's  law  and  equation  of  equilibrium 

=  ar'  +  r  (d<7r'/dr)  (27) 

can  be  used  to  express  az  in  terms  of  ar' 

az'  =  E£z'  +  2var '  +  vr(dcrr'/dr)  (28) 

Since  the  dilation  is  purely  elastic, 

du'/dr  +  u'/r  +  £z'  =  E-1  (l-2i>)  (or'  +oq'  +oz' ) 

Integrating  the  above  equation  after  substituting  from  Eqs.  (27)  and  (28), 
we  have 


ru'  =  E-1  (l-2v)  (l+t»)r2ar'  -  v£z'r2  +  c'  (29) 

where  c'  is  a  constant.  Now  the  strain  components  £0 ' ,  er ' ,  e0'e,  er'e  can 
all  be  expressed  in  terms  of  ez',  or',  and  d a' /dr. 

Since  €0'P  =  e0'  -  £0'e,  the  plastic  strain  components  are 

£0 1 P  =  -£r'P  =  -E-' (l-v2)r(dar'/dr)  +  c'/r*  (30) 

and  the  equivalent  plastic  strain  is 

-2  -2 

£ 1 P  =  --  £0'P  =  —  [c'/r^-E-1 (l-v«)rdar’/dr]  (31) 

/3  /3 

it  has  been  assumed  that  the  sign  of  d£0'P  is  the  same  throughout  the 
unloading  process  and  that  it  is  negative.  This  will  be  the  case  when  the 
internal  pressure  is  removed  during  the  unloading.  In  the  plastic  region,  a  $ 
r  ^  p' ,  we  have 

rdCTr'/dr  =  a0'  -  ar'  =  a'  (32) 


Substituting  Eq.  (32)  into  Eq.  (31),  we  can  determine  the  constant  c'  by  the 
condition  that  e'P  =  0  in  the  elastic  region,  p*  <  r  <  b,  and  obtain 


E_1  (1-u*  )  [(p'  2/r2  )<j0'| 

P' 


<7'1 


(33) 


With  the  aid  of  Eqs.  (3)  and  (4),  we  can  simplify  the  above  expression  to  the 
following  form: 


C  =  0i[g(fp')p,2/r2  -  g(C) ]  (34) 

where  g ( Cp • )  is  the  value  of  g(f)  given  by  Eq.  (4)  evaluated  at  r  =  p'  and 

/3  (1-m' ) 

*l'  "  (I-"*’  )/tm*  +  --  “'""I  (35) 


Equations  (32)  and  (33)  together  with  Eq.  (3)  suffice  to  determine  crr'  in  the 
plastic  region.  Integrating  Eq.  (32)  with  the  boundary  conditions  that  ar 1  is 
continuous  at  r  =  p'  and  ar'  =  p  at  r  =  a,  we  have 

°r'/°o  *  P/*o  ~  nCr * Ca)  (36) 
and 


ar ' /°o 


=  P/aQ  -  r(Cp’,Ca) 


P'  + 


(37) 


where 

1 

Hfr.Ca)  a  -  -  (^2,//51)(p,/p)2(Cr-Ca)g(Cp- )  +  ( 1~/J2 ’  )Q(Cr . Ca )  (38) 

fCa  j 

G(Cr»Ca)  =  Jf  -  (C^iJ-’gtOdC  (39) 

»r  2 

^2 '  =  m'jSi'/d-m' )  ,  g ( C )  is  given  by  Eq.  (4)  (40) 

and  Cr»  Ca  can  Pe  calculated  by  Eq.  (12).  For  some  special  cases,  e.g. 
Reference  10,  explicit  expressions  for  G({r,Ca)  can  be  obtained.  In  general, 
numerical  integration  is  needed. 
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The  stresses  in  the  elastic  zone  ( p '  $  r  <  b)  are 

or'/o0  !  (41) 

=  -  gKp'HifP'/r)*  -  (p'/b)*] 

Pe'/cr0  2  (42) 

The  other  expressions  in  the  entire  tube  (a  *  r  <  b)  are 

oz'  /a0  =  v(or'+<jQ')/o0  +  Eez'/cr0  (43) 

Eez'/a0  =  -  (u-2v) (P/a0)/(b*/a*-l )  (44) 


(E/o0)u'/r  =  (l-2v) (1+v) (Op'/Oo)  -  vEez’/a0  -  (l-v*) (p,/r)*g(Cp' )  (45) 

The  residual  stresses  and  the  residual  displacement  are  found  by  addition 
ar"  =  ar  +  ar'  ,  oq"  *  oq  +  eg '  ,  a2"  =  ctz  +  az' 

and  u"  =  u  +  u'  (46) 

NUMERICAL  RESULTS  AND  DISCUSSION 

The  numerical  results  for  three  closed-end  thick-walled  cylinders  with 
various  wall  ratios  have  been  obtained.  The  material  constants  in  all  cases 
were  v  =  0.3,  E/o0  *  200,  and  m  =  0.01.  The  Bauschinger  effect  factor  (f)  was 
determined  during  computations  based  on  the  experimental  curve  f(eP)  as  shown 
in  Figure  8  of  Reference  9.  The  slope  of  unloading  after  reverse  yielding  was 
estimated  to  be  m'  =  0.3  for  a  high  strength  steel.  The  computations  for  the 
stresses,  strains,  and  displacements  during  elastic-plastic  loading  and  after 
complete  unloading  were  made,  and  some  of  the  numerical  results  are  presented 
graphically  in  Figures  3  through  9. 

Figure  3  shows  the  relation  between  the  pressure  factor  and  the  dimen¬ 
sionless  elastic-plastic  boundary  for  a  tube  of  wall  ratio  two.  Also  shown  in 
the  figure  is  the  functional  relationship  between  the  residual  displacement  at 
the  bore  and  the  amount  of  overstrain.  The  new  results  for  the  residual 
displacement  shown  in  the  figure  are  a  little  smaller  than  the  well-known 
results  based  on  complete  elastic  unloading  or  the  isotropic  hardening  model. 


-p/a 


Figure  3.  Internal  pressure  and  residual  displacement  vs.  elastic-plastic 
interface  for  the  case  of  wall  ratio  two. 

A  similar  result  for  a  closed-end  tube  of  wall  ratio  three  is  shown  in  Figure 

4.  The  applied  pressure  to  reach  100  percent  overstrain  for  the  case  of  wall 

ratio  three  is  1.61  times  that  for  the  case  of  wall  ratio  two.  The  residual 

displacement  at  the  bore  has  been  increased  from  ua"/a  *  0.937  percent  for  the 

case  of  b/a  »  2  to  ua"/a  =  2.939  percent  for  the  case  of  b/a  *  3. 
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Figure  4.  Internal  pressure  and  residual  displacement  vs.  elastic-plastic 


interface  for  the  case  of  wall  ratio  three. 

Since  the  radial  and  axial  stresses  are  not  as  important  as  the  hoop 
stress,  only  the  numerical  results  for  the  latter  are  presented.  Figure  5 
shows  the  distribution  of  residual  hoop  stresses  in  a  tube  of  wall  ratio  two 
subjected  to  25,  50,  and  75  percent  overstrain.  The  elastic-plastic  locations 
and  the  pressure  factors  corresponding  to  three  cases  are  p/a  »  1.25,  1.50, 
1.75,  and  p/a0  *  0.5284,  0.6265,  0.6818.  Reverse  yielding  zones  occur  for  the 
second  and  third  cases,  but  not  for  the  first  case.  Figure  6  shows  the 

12 


distribution  of  residual  hoop  stresses  in  a  closed-end  tube  of  wall  ratio 
three  subjected  to  25,  50,  and  75  percent  overstrain.  Reverse  yielding  occurs 
for  all  cases  and  spreads  a  larger  zone  as  the  amount  of  overstrain  increases. 
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Figure  5.  Residual  stress  distribution  in  an  autof rettaged  tube  (b/a  =  2.0). 
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figure  6.  Residual  stress  distribution  in  an  autof rettaged  tube  (b  =  3a). 

It  is  of  interest  to  compare  the  new  results  for  residual  hoop  stresses 
with  those  based  on  the  other  theoretical  models.  Figure  7  shows  a  comparison 
of  residual  hoop  stresses  based  on  three  theoretical  models  for  an  auto- 
frettaged  tube  of  b/a  *  2,  p/a  =  1.75.  The  dotted  line  represents  the  well- 
known  solution  based  on  the  assumption  of  complete  elastic  unloading.  The 
same  solution  is  obtained  when  either  the  isotropic  or  the  kinematic  hardening 
model  is  used.  The  broken  line  represents  the  solution  based  on  a  model 
including  the  Bauschinger  effect  but  with  the  same  hardening  parameter  m'  =  m 


Figure  7.  A  comparison  of  residual  stress  based  on  three  material  models 
(b/a  =  2,  p/a  =  1.75). 

a  0.01.  This  solution  is  close  to  Parker's  solution  neglecting  hardening 
completely  with  m’  =  m  =  0  (ref  8).  As  can  be  seen  from  this  figure,  the  new 
result  represented  by  the  solid  line  is  quite  different  from  the  other  two 
solutions.  This  comparison  also  demonstrates  that  the  Bauschinger  effect  and 
strain-hardening  during  loading  and  unloading  can  have  a  significant  effect  or 
the  residual  stresses,  especially  near  the  bore.  For  a  closed-end  tube  of 
wall  ratio  three  subjected  to  50  and  100  percent  overstrain,  we  show  the 


comoarison  of  residual  hoop  stresses  based  on  three  models  in  Figures  8  and  9, 
respectively.  The  dotted  lines  represent  the  solutions  based  on  the  isotropic 
hardening  model  (ref  3).  Reverse  yielding  occurs  even  if  the  Bauschinger 
ef*e ct  is  neglected.  The  solid  and  broken  lines  represent  the  new  model  and 
the  model  with  m'  *  m,  respectively.  As  can  be  seen  from  these  figures,  the 
numerical  results  of  the  residual  hoop  stresses  based  on  three  models  are 
quite  different  in  the  reverse  yielding  zone.  The  new  result  is  considered 
more  accurate  for  our  problem  because  it  is  based  on  a  more  general  theoreti¬ 
cal  model  considering  the  Bauschinger  and  hardening  effects  during  loading  and 
unloading. 
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Figure  8.  A  comparison  of  residual  stress  based  on  three  material  models 
(b/a  =  3,  p/a  =  2.0). 
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Figure  9 
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A  comparison  of  residual  stress  based  on  three  material  models 
(b/a  =  3,  p/a  =  3.0). 
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